QCD equation of state in a virial expansion 
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We describe recent three-flavor QCD lattice data for the pressure, speed of sound and interaction 
measure at nonzero temperature and vanishing chemical potential within a virial expansion. For 
the deconfined phase we use a phenomenological model which includes non-perturbative effects from 
dimension two gluon condensates that reproduce the free energy of quenched QCD very well. The 
hadronic phase is parameterized by a generalized resonance-gas model. Furthermore, we extend 
this approach to finite quark densities introducing an explicit /^-dependence of the interaction. We 
calculate pressure, quark-number density, entropy and energy density and compare to results of 
lattice calculations. We, additionally, investigate the structure of the phase diagram by calculating 
the isobaric and isentropic lines as well as the critical endpoint in the (T, /iq)-plane. 
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I. INTRODUCTION 



The exploration of the phase structure of quan- 
tum chromodynamics (QCD) for the whole density- 
temperature plane is a challenging task for the standard 
model. Lattice Monte Carlo simulations have revealed 
several exciting results over the past decade [l|, [3, 0, H, 
H, Sill- At finite densities, modehng of QCD has added 
substantially to our understan ding of its rich phase struc- 
ture I, M El, M M M M, MM, MMMM- in 

the last years heavy-ion collision experiments focus on 
the investigation the quark-gluon plasma (QGP). Early 
concepts of the QGP were guided by the assumption of 
a weakly interacting systems of partons because the en- 
tropy density s and the energy density e calculated from 
lattice QCD were close to the Stefan Boltzmann (SB) 
limit for a relativistic noninteracting system. However, 
recent observations at the Relativistic Heavy Ion Col- 
lider (RHIC) indicated that the QGP created in ultrarel- 
ativistic Au + Au collision was interact ing more strongly 
than hadronic matter [H [H, [H, HI, Hi 113, HI . Be- 
cause of these experimental facts the noticeable deviation 
(w 10 — 15%) became more important. 

In the hadronic phase of QCD the partition function of 
a hadron resonance-gas yields a satisfactory description 
of lattice results on bulk thermodynamic observables in 
the low temperature region for = (29| and in the 
general case at finite densities, i.e. ^ [s^l- In this 
context the investigation of the interaction between the 
partons in the deconfined phase plays a crucial role to 
understand the properties of the QGP. The dynamical 
quasiparticle picture is one important way to relax the 
approximation of the QGP as an ideal gas [itI. ITsl [Sll. [3^ 
of massive but noninteracting degrees of freedom. Alter- 
natively, the phase structure can be investigated in the 
finite width model (FWM) by analizing its isobaric par- 
tition given by a generalization of statistical ensembles 
with fluctuating extensive quantities (ssl . Is^ ]. Another 
possibility is to use a generalized formalism to calcu- 
late the corrections to a single-particle partition function 
starting from an interaction potential. 



In this work we use the latter formalism, which is 
a generalized version of the classical virial expansion, 
to calculate the partition function with an interaction 
potential extracted from lattice calculations. 

The paper is organized as follows: In Section |TT] we 
introduce the virial expansion for vanishing quark chem- 
ical potential /iq. In Section UIIl we motivate our model 
for the hadronic as well as for the deconfined phase and 
in Section IIVI we show the results of the virial expan- 
sion within our model. We extend in Section |V] the virial 
expansion to finite densities and again compare our re- 
sults with lattice calculations. Furthermore, we explicitly 
evaluate the confinement-deconfinement phase transition 
and the isobaric and isentropic lines. The conclusions in 
Section IVTl finalize this work. 



II. VIRIAL EXPANSION AT /i, = 

The main quantity to establish thermodynamical prop- 
erties of QCD is the partition function Z{T,V). In this 
Section we restrict ourselves to the case of vanishing 
chemical potential /iq, i.e. vanishing baryon number. The 
partition function then is obtained as 



Z{T, V) = Tr [e-'^^j 



(1) 



where H is the Hamiltonian of the system and f3 = 1 /T is 
the inverse temperature. For particles with spin-isospin 
degeneracy factor d the logarithm of the partition func- 
tion reads 



-f3H{v,p] 



(2) 



where 77 = —1 for bosons and 77 = 1 for fermions. 

In the case of an interaction free system, the Hamilto- 
nian H{r, p) is given by the kinetic energy Ek and there- 
fore depends on p = |p| only. For massless particles one 
has (with h — c = 1) H{p) ^ p and can calculate the 
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partition function for bosons as well as for fermions, 

77r2 



The pressure reads 



InZp = dpV 



720/33 F 



(0) 



(4) 



Accordingly, the boson and the fermion partition func- 
tions are connected by 



dBln4°^ = ldF\nZ^\ 



(5) 



For non-interacting particles with mass m the Haniilto- 
nian is given by H{p) — y'-p^ + rn^ and the partition 
function can be expressed in terms of the Bessel function 

K2. 

From the partition function all thermodynamic quanti- 
ties can be derived: pressure, energy density, entropy, in- 
teraction measure and velocity of sound. Using the free 
massless partition function this leads to the well-known 
Stefan Boltzmann limit. 

To relax the dependence on the size of the system we 
divide the partition function by the volume V and all 
calculated quantities become densities. In general we are 
not able to calculate exactly the partition function except 
by numerical diagonalization of the many-body Hamilto- 
nian. To incorporate the effect of interactions we start 
with the classical virial expansion. To this aim we con- 
sider the A^-body Hamiltonian defined by 



^fi(i"i,Pi) 
-ff2(ri,r2,Pi,P2) 
i?3(ri,r2,r3,pi,p2,P3) 



= E^{p,) (6) 

= E^{pi) + E^{p2) 

+Ei,{p3) + W123 (8) 



In classical physics the Bose as well the Fermi distribu- 
tions in Eq. ^ merge to the Boltzmann distribution and 
the logarithm in the integrand drops. Then, we can sepa- 
rate the kinetic and the interaction term in the partition 
function. To recover the quantum statistics we perform 
a minimal substitution by replacing the integral over the 
momentum by the integral for the free case, which leads 
to the exact Stefan Boltzmann limit. The definition 



C = lnZ("), 



(9) 



where we suppressed for simplicity the indices for fermion 
and bosons, leads to an expansion of the partition func- 
tion in powers of ^, i.e. 



(10) 



u=0 



We consider only the first two terms of this expansion: 
then the coefficients are given by 



bi 



1 



(11) 
(12) 



p = rinZ = TV4c 

7/1 



(13) 



and can be directly expressed in terms of the particle 
density p. We obtain for the first terms of the expansion 



(14) 



The other quantities can be calculated from the pressure 
or the partition function using the thermodynamic rela- 
tions: 



dlnZ 

^ ~ w~ 

dP 

o dP dP/e P 

Cs = -J- ^^—T— + —■ 

de de e 



(15) 
(16) 
(17) 



III. THE MODEL 

The application of the virial expansion formalism to 
investigate the properties of the hadronic and partonic 
matter requires the calculation of the single-particle par- 
tition function for the degrees of freedom and of the cor- 
rections given by the residual interaction between the 
'constituents'. For the confined phase the degrees of free- 
dom are the hadrons while in the QGP phase they are 
the quarks/antiquarks as well as the gluons. 

The analysis of thermodynamic properties of the 
hadronic phase has demonstrated that a good descrip- 
tion is achieved using a phenomenological resonance-gas 
model taking into account all mesonic and baryonic res- 
onances with masses up to 1.8 GeV and 2.0 GeV, respec- 
tively. This amounts to 1026 resonances that truly are 
experimental parameters of the model [S?]. In this case 
the residual interaction may be neglected. However, to 
reduce the number of the experimental parameters for 
the description of the low temperature region of QCD we 
explicitly write the partition function of the resonance- 
gas in the Boltzmann limit as 



35| 



inz(y,r) = -^P^^^)K^ [-y) ' 



(18) 



where the sum begins with the stable ground state mg 
and includes the resonances m^, i = 1, 2, • • • with weights 
p{mi) relative to ttiq. To specify p{mi) we employ the 
Hagedorn hypothesis [3^ [3^ of an exponential growth 
of the number of hadronic resonances with mass. The 
formula for the asymptotic dependence of the density of 
hadronic states on mass m is given by 



p{m) = /(m)exp{m/rH}, 



(19) 
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where f{m) denotes a slowly varying function and Th is 
the Hagedorn temperature estimated to be in the range 
of 180 - 210 MeV. We choose the form 

f{m)^am-°', a 3/2 (20) 

defined for mo > 140 MeV taking into account that 
the pions are the hadronic states with the lowest mass. 
We now replace the sum in the resonance-gas partition 
function ()18|) by an integral and insert the exponentially 
growing mass spectrum (fT9|) with the choice of /(m) ([20]) . 

aVT f°° /m\ 
In Z{V, T) = / dm V^exp {m/THji-Ta (^) • 

(21) 

In the limit of low temperatures, i.e. z — m/T 3> 1, the 
Bessel function can be approximated by, 

^2(^) = y^exp{-z}, (22) 

and the integral in the resonance-gas partition function 
([2T]l can be solved analytically. This gives 

^ ' ' \2tt) b ' T Th 

(23) 

and leads to a good description of the thermodynamical 
variables. 

To perform a realistic investigation of the properties 
of the QGP we have to determine a realistic effective in- 
teraction between the partons in the deconfined phase. 
The starting point for this aim is provided by lattice 
calculations for the quark-antiquark free-energy. There 
are many parametrizations of this quantity and several 
discussions about the role as effecti ve p otentials (or al- 
ternatively as the internal energy) [33]. We adopt the 
parametrization of the free energy given by Megias et 
al. in Ref. [s^. This phenomenological model includes 
non-perturbative effects from dimension two gluon con- 
densates. This one-giuon exchange model, introduced in 
Ref. [40| for the description of the Polyakov loop, allows 
for a good parametrization of the singlet heavy quark- 
antiquark free energy lattice data in gluon-dynamics. 
The connection between dimension two gluon conden- 
sates and QGP-properties has been confirmed by recent 
fits of the interaction measure in gluon-dynamics in terms 
of the condensate [13, 13 ■ 

The expression for the free energy in this model is given 

by 

F,{r,T)^Vi{r,T) + F,{^,T), (24) 

where ^1(00, T) is the singlet free energy at infinite sepa- 
ration and Vi is the effective potential for the deconfined 
phase. We use 

V,{r,T) = -(-- + -^)e'^"^^>, (25) 
^ ' ' \12r 2NcT J ' ^ ^ 



where C2 is the non-perturbative dimension two conden- 
sate and M{T) a Debye mass estimated as 

M{T) = V^VS+We gT = gT. (26) 

The parameters have been determined in Ref. 39] from 
a fit to the lattice data of the Polyakov loop in gluon- 
dynamics, i.e. N[ = 0, 

C2 = (0.9GeV)2 (27) 
gpL = 1.26. (28) 

In our calculations we will consider light quarks instead 
of heavy quarks (with infinite mass) and hence expect 
deviations from the previous values (|27p and (pS)) . Ac- 
cordingly, we consider ^ as a free parameter that we may 
fix in comparison with A'f = 3 lattice calculations 0] for 
the pressure (calculated for Nt = 6). These lattice cal- 
culations imply a critical temperature Tc ~ 196 MeV for 
the deconfinement phase transition. 

IV. RESULTS 

Using the one-gluon exchange potential we can directly 
calculate the pressure from Eq. (fT5)) for the QGP. To 
calculate the second virial coefficient 62 we linearize the 
exponential function in Eq. (|12p . For the hadronic phase 
we calculate the pressure from the partition function us- 
ing Eq. By comparing the lattice results of the 
pressure P/T^ as a function of the temperature for the 
deconfined phase (given by T > Tc = 196 MeV) with our 
results we find a good agreement with g in the range 
1.25 < g < 1.35. Note that these values for g are 
comparable with the result from the fit of the Polyakov 
loop, i.e. gpL — 1.26. Our fit in the deconfined region 
yields g — 1.30 and we can fix the value of the Hagedorn 
temperature Th and the constant a to obtain a contin- 
uous function for all temperatures T > 0. This leads to 
Th = 205 MeV and a = 9.17 MeV^^^ 

In Fig. [1] we show the lattice results of the scaled pres- 
sure P/T^ as a function of the temperature expressed in 
units of the critical temperature Tc = 196 MeV and the 
pressure from the virial expansion with the parameters 
fixed above. 

Since all thermodynamic quantities are connected by 
known relations (cf. (|14p - pT|) ) we expect that the lattice 
data for all quantities are compatible with the results 
from our model, too. To confirm this we show in Fig. [2] 
the entropy density s divided by as a function of the 
temperature expressed in units of the critical tempera- 
ture using Eq. (|16p . 

Another interesting quantity is the square of the veloc- 
ity of sound defined by Eq. pT)) . In the high temper- 
ature limit as well as in the transition region, where the 
derivative d{P/e)/de vanishes, Cg is given by p/e. The re- 
sults of our model for this quantity (as a function of the 
fourth root of the energy density) are shown in Fig. [3] for 
Tc = 196 MeV (solid) in comparison with the lattice data 
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(blue squares). They show (as expected) a good agree- 
ment with the lattice data for e'^/^ / [GeV /im^f/^ > 1.3 
which is equivalent to T > Tc- The dashed line is the 
parametrization in the high temperature region with the 
simple Ansatz 



1 



-A ) 

Beim^/GeYj 



(29) 



that provides a good fit to the data in the interval 1.3 < 
ei/V(GeV/fm3)i/4 < 5 with C = 0.964(5), A = 1.16(6) 
and B = 0.26(3) 3. The current resolution and accuracy 
of lattice calculations for the hadronic region clearly are 
not yet sufficient to allow for a detailed comparison with 
any model calculation. The dotted line in Fig.[2]indicates 
the exact square of the sound velocity. The difference 
between this quantity and the ratio P/e indicates the 
deviation from a linear dependence between pressure and 
energy density valid at high temperature. This linear 
behavior is confirmed in Fig. [3] for high energy densities, 
where Cg as well as P/e show the Stefan Boltzmann limit 
of 1/3. 

In the framework of the one-gluon exchange model the 
interaction measure W = e — 3P plays an important role 
and was recently investigated in detail in Refs. [13, 13 ■ 
From the fit of W in the deconfined region in gluon- 
dynamics the value of the dimension two gluon conden- 
sate has been calculated in an independent way without 
considering the Polyakov loop. However, the results of 
the two approaches show only little deviation when us- 
ing 1.26 < g < 1.45. The calculation of the interaction 
measure using the virial expansion starting with a model 
motivated by the fit to the Polyakov loop may directly 
indicate a dynamical connection between the two meth- 
ods. For this reason we show in Fig. 2] the interaction 
measure divided by the fourth power of the temperature 
W/T"^ as a function of the temperature expressed in units 




FIG. 1: (Color online) Pressure P as a function of the 
temperature divided by for g = 1.30, Th = 205 MeV 
and a = 9.17 MeV^^^ in comparison to the lattice data 
from Ref. M. 
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FIG. 2: (Color online) Entropy density s as a function of 
the temperature divided by for g = 1.30, Th = 205 
MeV and a = 9.17 MeV^''^ in comparison to the lattice 
data from Ref. 
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FIG. 3: (Color online) Ratio P/e as a function of the 
fourth root of the energy density for the critical temper- 
ature Tc = 196 MeV in comparison to the lattice data 
from Ref %. 



of the critical temperature Tc = 196 MeV, where the ex- 
pected agreement with the lattice data (blue squares) is 
observed. 



V. FINITE QUARK CHEMICAL POTENTIAL 

The extension of the virial expansion to non- vanishing 
quark chemical potentials [iq requires an explicit depen- 
dence of the interaction on the quark chemical potential. 
Since the interaction is primarily determinated by the 
Debye mass M, we may generalize the interaction for 
non-vanishing [iq by introducing an explicit dependence 
for the Debye mass only, i.e. M(T) — > M{T,fiq). 
We parameterize this dependence by 



M(T,/.,) = M(")(T)P(^) 



(30) 
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FIG. 4: (Color online) The interaction measure W as a 
function of the temperature divided by for g — 1.30, 
Th = 205 MeV and a = 9.17 MeV^''^ in comparison to 
the lattice data from Ref. [gl. 



where Af(")(T) = M{T, Hq = 0) denotes the Debye mass 
at vanishing chemical potential and F describes the mod- 
ification at finite density. This assumption is equivalent 
to introduce a /iq-dependence of the coupling g (in our 
model) , 
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9 



(31) 



Furthermore, perturbation theory suggests that the lead- 
ing order contribution to the correction is given by (4^ 



3iVf 



(2iVe -I- Nt) 



(32) 



For our calculations we generalize this expression as fol- 
lows. 



= A 



with 



A 



V T 



(33) 



(34) 



where the value of the coefficient tt^^ is consistent for 
A^f = 3 with Ref. fZsj. For the hadronic phase we use an 
effective modification of the constant a given formally in 
Eg. ipTj) . Choosing in this case the perturbative correc- 
tion we obtain 



(35) 



where /is = 3/ig indicates the baryonic chemical poten- 
tial. 

At nonzero chemical potential, typical quantities of in- 
terest that become accessible from lattice calculations are 
the (scaled) pressure difference defined as 



AP(r,/i,) _ p{T,^lq)-PiT,^lq^o) 



rp4 



rp4 



(36) 



and the scaled quark number density given by 



T3 



(37) 



Since for [iq^Q direct Monte Carlo calculations are not 
applicable, higher order derivatives of the pressure with 
respect to [iqjT at muq = are calculated and then the 
pressure is estimated using a Taylor expansion, 



p(T,A.,) _ inz(r,/i,) _^ mr^r 



(38) 



where the expansion coefficients are given in terms of 
derivatives of lnZ(T, ^q), i.e. 



c„(T) 



1 9" In Z 



(39) 



The series is even in iiq/T which reflects the invari- 
ance of Z{T, ^q) under exchange of particles and antipar- 
ticles. The coefficients C2 and C4 have been calculated for 
two-flavor QCD in Ref. [1]; an extension to the sixth 
order is given in Ref. Q. To evaluate the pressure and 
the other thermodynamics quantities within the virial ex- 
pansion we can formally use the expansion given in Eq. 
psp . where the expansion parameter C now depends on 
the quark chemical potential, i.e. C, — > Cit^q)- Using the 
Taylor expansion of the pressure in the Stefan Boltzmann 
limit we can explicitly give this /i^-dependence of (. In 
this case the coefficients are well known as 



„SB 



Nf 

T' 



SB _ ^2 



i2ny 



c^^ = for n>6. 



(40) 

Now we can rewrite the generalized expansion parameter 
as 



(41) 
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FIG. 5: (Color online) Scaled pressure difference as a 
function of the temperature for different values of i-iq/T 
compared to lattice calculations taken from Ref. (Jj. 
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FIG. 6: (Color online) Scaled pressure difference as a FIG. 8: (Color online) Coefficient C4{T) of the Taylor 

function of the temperature for different values of fiq com- expansion for the pressure as a function of T/Tc compared 

pared to lattice calculations taken from Ref. to lattice calculations taken from Ref. 0]. 



The results for the scaled pressure difference as a function 
of the temperature (expressed in units of the critical tem- 
perature Tc = 196 MeV) are shown in Fig. O for various 
values of ^J.q/T. The data points are the lattice results 
for N{ = 2 from Ref. [3] scaled to iVf = 3. Lattice results 
for three-flavor QCD are not yet available a nd p artial 
data for A^f = 2 -I- 1 are preliminary only [H, HI] . The 
agreement between our results and the lattice simulations 
is quite satisfactory. In the region around Tc there is a 
deviation from the scaled lattice points. 

This is better visible in an analogue representation of 
the results in Fig. [51 where the scaled pressure difference 
is shown as a function of the temperature expressed in 
units of the critical temperature = 196 MeV for var- 
ious values of the quark chemical potential fig. We note 
an overestimation of the pressure difference around the 
critical temperature. This can be understood by consid- 
ering the possible errors of the scaling procedure as well 



as by the errors given by the (truncated) Taylor expan- 
sion itself. 

We can now calculate the coefficients C2(T) and C4{T) 
of the Taylor expansion for the pressure by the least- 
square method. The comparison with the scaled lattice 
results from Ref. |4] , shown in Fig. [H is quite satisfac- 
tory. Analogue comparisons with quasiparticle models 
are given in Refs. [46l.l47t. Although these coefficients al- 
low to directly evaluate the scaled quark number density, 
we prefer to use the definition given in Eq. ([37|) because 
the Taylor sum is only an expansion up to the fourth 
order and its validity for strongly interacting system is 
not so clear (cf. Ref. [50]). Our results for Uq/T^ as 
a function of the temperature (in units of the critical 
temperature Tc = 196 MeV) are shown in Fig. [5] for var- 
ious values of the quark chemical potential /ig. Also in 
this case, the comparison between our virial calculation 
and the corresponding (scaled) lattice data shows a slight 
overestimation close to Tr. 
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FIG. 7: (Color online) Coefficient C2(r) of the Taylor 
expansion for the pressure as function of T/T compared 
to lattice calculations taken from Ref. [J]. 



FIG. 9: (Color online) Scaled quark number density as 
a function of the temperature for different values of 
compared to lattice calculations taken from Ref. 
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Using the results for the pressure and the entropy den- 
sity (from Eg. pB)) ) we can directly investigate the struc- 
ture of the QCD phase diagram in the fiq — T plane. 
Starting from the point at zero chemical potential, i.e. 
(0, Tc), we calculate the isobaric {P — const.) as well as 
isentropic (s = const.) lines in the phase diagram which 
are shown in Fig \W\ by the solid and dot-dashed lines, 
respectively. The resulting isentropic temperature at the 
phase boundary is practically constant up Hq « 200 MeV, 
whereas the isobaric temperature at the phase bound- 
ary (solid line) drops with increasing chemical potential. 
This behavior quantitatively agrees with the critical line 
(dashed) from lattice calculations within the reweighting 
method from Ref . [5l| . The observation that the critical 
phase transition is a line at constant pressure is also con- 
firmed by lattice simulations in Ref. [51]. The dotted line 
in Fig. 10 shows the region fi < T, where the lattice re- 
sults are expected to be valid. For comparison we display 
also the hadron chemical freeze-out line (dash-dash-dot) 
using the parametrization given in Ref. |52l |. Note that 
the hadron freeze-out - determined experimentally from 
hadron ratios - clearly occurs in the hadronic phase for 
all fiq with a sizeable gap in temperature relative to the 
phase transition line. 

Furthermore, a detailed investigation of the isobaric 
line leads to a first estimate for the critical endpoint 
{iie,Te), where the crossover in the low density re- 
gion changes to a first order phase transition. Following 
Ref. the curvature of the crossover - separating the 
QGP and the hadronic phase - is given by 



200 - 



T/T, = 1 - Cf?q/T^, 



(42) 



To evaluate the critical endpoint we perform a corre- 
sponding parametrization of our calculated isobaric line. 
This only works if ^ 1j i-^- ^ maximun of the chem- 
ical potential exists which satisfies this condition; the 
latter value defines the chemical potential of the end- 
point. From our analysis we obtain /i^ = 117±5, which 
is in agreement with the results of the lattice calculation 
^{Latuce) ^ j^20 ± 10 from Ref. The corresponding 

critical point (fiEjTs) — (117,190) MeV is shown as a 
star in FiglTOl 

We recall that the main aim of this paper is to evaluate 
the QCD equation of state in terms of physically mea- 
surable quantities. Following Ref. [4| we calculate the 
equation of state A(p/r^) vs (n^/T^)^ for different tem- 
peratures. The results are shown in Fig. [TTl where the 
yellow region is delimited by the low-T SB (upper bound- 
ary) and by the high-T SB limit (lower boundary). For 
temperatures near- but below - Tc the equation of state 
shows a form resembling the low-T SB limit 
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FIG. 10: (Color online) Overview of the QCD phase dia- 
gram with isobaric (solid) and isentropic lines (dash-dot) 
in comparison to the phase transition line (dashed) esti- 
mated in Ref. [sJl. Also shown is the hadron chemical 
freeze-out line (dot-dash-dashed) in the parametrization 
given in Ref. j52|. The dotted line (T = /iq) demonstrates 
the expected range of validity for the IQCD calculations. 
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FIG. 11: (Color online) Equation of state /^{p/T'^) vs 
/r^)2 for various temperatures. The continuum Ste- 
fan Boltzmann forms for low and high T (given by (|43| 
and (|44p ) are also shown as functions of T/Tc. 



form of the high-T SB limit given by 



(44) 



The quadratic coefficient depends smoothly on the tem- 
perature and therefore shows a sizeable deviation from 
the exact Stefan Boltzmann limit also for T = Tn. 



(43) 



VI. CONCLUSIONS 



As T increases though the critical temperature, the EoS 
changes from this functional dependence to the stiffer 



We have performed an investigation of QCD thermo- 
dynamic quantities within a virial expansion approach. 
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While hadronic matter is described by a generalized 
resonance-gas model, we fix the interaction in the QGP 
phase by lattice calculations for the heavy quark free en- 
ergy. We adopt a phenomenological model which includes 
non-perturbative effects from dimension two gluon con- 
densates that describe the free energy of quenched QCD 
very well at vanishing chemical potential. Explicit ex- 
pressions for the partition function and the pressure in 
the deconfined phase are derived. Our results for fiq — 
and N{ = 3 for pressure, entropy density, speed of sound 
and interaction measure at nonzero temperature com- 
pare well with three-flavor QCD lattice calculations with 
almost physical masses from Ref. Q . We find that a cou- 
pling parameter g = 1.30 provides a good agreement for 
all quantities. Note that this choice is close to the value 
determined in Ref. [39] from a fit to lattice data for the 
Polyakov loop in gluon-dynamics, i.e. ^pl = 1.26. The 
difference in the coupling parameters most likely is due 
to the different quark masses employed. 

Furthermore, we have extended this approach to finite 
densities introducing an explicit /x-dependence of the in- 
teraction (via the Debye mass following perturbation the- 
ory). This allows to calculate the thermodynamic quan- 
tities of interest such as pressure and quark-number den- 
sity etc. also at finite /z,; the results again compare well 
with lattice calculations scaled from iVf = 2 to iVf = 3 
with minimal deviations at the critical temperature Tc- 
These differences are most visible in the pressure, become 
smaller for the scaled quark number density and lead to 
the open question about the validity of the Taylor expan- 
sion and its thermodynamic consistency. 

Additionally, we have investigated the structure of the 
phase diagram also in the region inaccessible to lattice 



calculations, i.e. /Lt, > T. We have calculated the isen- 
tropic and specially the isobaric line, which agrees re- 
markably well with the cross over transition evaluated 
by lattice calculations with the rcwcighting method. The 
direct knowledge of the partition function allows to cal- 
culate the equation of state of the quark gluon plasma 
which also reproduces very well the predictions of the 
lattice calculations. The EoS clearly shows two decisive 
properties of the QCD phase diagram: i) There is an 
unequivocal transition from the hadronic matter to de- 
confined matter as is evident from the change of the func- 
tional form of the equation state from the typical Stefan 
Boltzmann behavior for low to high temperatures, ii) 
The QGP is a strongly correlated system and sensitively 
differs from the ideal gas limit. Also at T = 2Tc there is 
a remarkable deviation of the EoS of the QGP from the 
Stefan Boltzmann limit which is a clear signal of surviv- 
ing correlations in the deconfined phase. 

Our equation of state has been dynamically calculated 
with a particular form of the interaction. Using the same 
interaction in molecular dynamical calculations we may 
systematically investigate relevant properties of the QGP 
such as the pair correlation function g{r), the plasma pa- 
rameter F, the structure factor S{q) or the shear viscos- 
ity rj. Further interesting applications in this context (cf. 
Refs. [13, [5^) are relativistic molecular dynamic simula- 
tions, where the phenomenological quark-quark interac- 
tion used here can be implemented as well as the QCD 
equation of state. This will allow to study partonic sys- 
tems also out of equilibrium. 

Acknowledgements: Work supported by the Deutsche 
Forschungsgemeinschaft (DFG). 
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